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Measurement of i^*^ production 
and the charm contribution to F2 
in deep inelastic scattering at HERA 



ZEUS Collaboration 



Abstract 

The production of L>*±(2010) mesons in deep inelastic scattering has been measured 
in the ZEUS detector at HERA using an integrated luminosity of 37 pb~^. The decay 
channels D*^ L'°7r+ (+c.c.), with —>■ K^ir^ or i^^7r"7r+7r+, have been 

used to identify the D mesons. The e'^p cross section for inclusive D*^ production with 
1< Q2 < 600GeV2 and 0.02 < y < 0.7 is 8.31±0.31(stat.)+o;5o(syst.) nb in the kinematic 
region 1.5 < pt{D*^) < 15GeV and \ r]{D*^)\ < 1.5. Differential cross sections are 
consistent with a next-to-leading-order perturbative-QCD calculation when using charm- 
fragmentation models which take into account the interaction of the charm quark with the 
proton remnant. The observed cross section is extrapolated to the full kinematic region 
in pt{D*^) and r){D*^) in order to determine the charm contribution, F2'^(x, Q^), to the 
proton structure function. The ratio ¥2"^ /F2 rises from ~10% at ~ 1.8 GeV^ to ~30% 
at ~130 GeV^ for x values in the range 10"'' to 10^^. 
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1 Introduction 



The first HERA measurements of tlie cliarm contribution, F2'^, to tlie proton structure function 
F2 were reported by tlie HI and ZEUS collaborations from the analyses of their 1994 deep 
inelastic scattering (DIS) data sets ||]. These early results, which were statistically limited, 
revealed a steep rise of as Bjorken-x decreases. At the lowest accessible x values, it was 
found that around 25% of DIS events contained open charm, in contrast to the EMC fixed 
target measurements in the high-x region where the charm contribution is small. Given the 
large charm content, the correct theoretical treatment of charm for F2 analyses in the HERA 
regime has become essential. More detailed measurements of charm production will aid such 
analyses. 

The early results suggested that the production dynamics of charmed mesons in ep 
collisions are dominated by the boson-gluon-fusion (BGF) mechanism shown in Fig. |I[ In this 
case, the reactions e^p e^D*^X are sensitive to the gluon distribution in the proton 
The measurement of charm production can also provide tests of perturbative QCD (pQCD), 
in particular, tests of the hard scattering factorization theorem, which states that the same, 
universal, gluon distribution should contribute to both F2 and F2^. In addition, the presence 
of two large scales, namely, the virtuality of the exchanged boson (Q^) and the square of the 
charm-quark mass {ml), provides a testing ground for resummation techniques. 

This paper reports a measurement of Z}*^(2010) production using the 1996 and 1997 data 
sets collected with the ZEUS detector, corresponding to an integrated luminosity of 37 pb~^. 
During this period, HERA collided E^. = 27.5 GeV positrons with Ep = 820 GeV protons, 
yielding a center-of-mass energy, ^/s, of 300 GeV. This larger data sample, together with 
some improvements to the ZEUS detector, allows an extension of the kinematic range to both 
larger and smaller Q^. The D*^ production is investigated in the decay channel D*+(2010) — >■ 
D'^(1864)7r+ (+ c.c), where 7r+ refers to a slow 7r+ f^, followed by the decay channels 
—>■ K^TT^ (+ c.c.) or —>■ K^ti^ti^ti^ (+ c.c). These two final states will be referred to 
as the K27T and KAn channels, respectively. 

Following a brief discussion of three pQCD calculations for D*^ production in DIS, a descrip- 
tion of the experiment and details of the data analysis are provided. Next, the cross sections for 
D*^ are presented and are compared to the QCD predictions. A possible cause of the observed 
discrepancies between the data and the QCD results is discussed. Finally, the D*^ results are 
used to obtain the charm contribution to F2. 

2 Charm quark production models 

General agreement was found between the earlier HERA D*^ data and the results from a Monte 
Carlo program HVQDIS for the production of heavy quarks and their subsequent fragmen- 
tation to hadrons. This program is based on next-to-leading-order (NLO) calculations 0] of 
order in the coefficient functions in the so-called fixed-fiavor-number scheme (FENS). In 
this scheme, the number of active quark flavors is fixed, independent of Q^. Only light quarks 
(m, d, s) are included in the initial-state proton and charm quarks are produced exclusively 
by BGF with NLO processes @. The presence of the two large scales, and m^, can spoil 
the convergence of the perturbative series because the neglected terms of orders higher than 
ttg contain log(Q^/m^) factors that can become large. Therefore the results of HVQDIS are 
expected to be most accurate at ~ and to become less rehable when ^ m^. 

In this high region, the c quark can be treated as massless, as implemented in the zero- 
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mass, variable-flavor-number scheme (ZM-VFNS). In this scheme, the resummation of large 
logarithms of Q"^ /ml |T0| results in a charm density which is added as a fourth flavor and 
which is then evolved in the same way as the light quark densities. At intermediate values, 
the two schemes need to be merged. One way in which this is done is described by ACOT 0] 



and by Collins |]TT|. An alternative matching method has been proposed by MRST |T2[ . 
A third method for modelling D*^ production has recently been suggested by BKL |13 



This tree-level pQCD calculation, applicable for > rric, considers the hadronization of the 
(cg)-state into a D*, in contrast to hadronizing an isolated c-quark. The D* is created from both 
color singlet and color octet configurations of the light and heavy quarks. Results from e'^e~ 
annihilation imply that the octet contribution is small. However, the singlet contribution alone 
underestimates the ZEUS data on the photoproduction of charm |14]. This calculation has 



been extended to DIS charm production and is compared to the Z)*^ data reported here. 

3 Experimental setup 



ZEUS is a multipurpose detector which has been described in detail elsewhere |T^. The key 



component for this analysis is the central tracking detector (CTD) which operates in a magnetic 
field of 1.43 T provided by a thin superconducting solenoid. The CTD ||T6| is a drift chamber 
consisting of 72 cylindrical layers, arranged in 9 superlayers covering the polar angleQ region 
15° < 6 < 164°. The transverse momentum resolution for full-length tracks is (t{pt)/pt = 
0.0058 0.0065 0.0014/pt {pt in GeV). The CTD was also used to establish an interaction 
vertex for each event. 

The uranium-scintillator sampling calorimeter (CAL) surrounds the solenoid. The CAL is 
hermetic and consists of 5918 cells, each read out by two photomultiplier tubes. The CAL 
contains three parts, the forward (FCAL), barrel (BCAL) and rear (RCAL), with longitudinal 
segmentation into electromagnetic and hadronic sections. The energy resolutions, as measured 
in test beams, are a/E = 0.18/ \J E{GeM) and 0.35/-y/-E'(GeV) for electrons and hadrons, re- 
spectively [jl7 . 



The position of positrons scattered close to the positron beam direction is determined by 
a scintillator strip detector (SRTD) fl^. The luminosity was measured from the rate of the 
bremsstrahlung process, e^p —>■ e'^'jp, where the photon is measured by a lead/scintillator 



calorimeter 19 located at Z = — 107 m in the HERA tunnel. 



4 Kinematics and reconstruction of variables 

The reaction e~^{k) + p{P) e^{k') + X at fixed squared center-of-mass energy, s = {k + P)^, 
is described in terms of = —q^ = —{k — k'Y and Bjorken-x = Q'^/{2P ■ q). The fractional 
energy transferred to the proton in its rest frame is y = Q'^/{sx). The virtual photon (7*)- 
proton center-of-mass energy W, given by W'^ = {q + -P)^, is also used, see Fig. |]. 

In neutral current (NC) e~^p DIS, both the final-state positron, with energy E'^ and angle 
6'g, and the hadronic system (with a characteristic angle 'jh, which, in the simple quark-parton 
model, is the polar angle of the struck quark) can be measured. The scattered positron was 



identified using an algorithm based on a neural network pO| . CAL cells were combined to form 



^ The ZEUS coordinate system is defined as right-handed with the Z axis pointing in the forward (proton 
beam) direction and the X axis horizontal pointing towards the center of HERA. The origin is at the nominal 
ep interaction point. The polar angle 9 is defined with respect to the positive Z direction. 
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clusters and combinations of these clusters and CTD tracks were used to reconstruct energy- 
flow objects (EFO's) [0, For perfect detector resolution and acceptance, the quantity 
S = — Pz,i) is equal to 2Ee (55 GeV). Here, Ei and p^^i are the energy and longitudinal 

component of the momentum assigned to the i-th EFO. The sum runs over all EFO's including 
those assigned to the scattered positron. 

In the K271 analysis, was reconstructed from the scattered positron {Ql) with the electron 
method 03] and y with the so-called S method [El] 



= (1) 

where 6had is calculated in the same way as 6 but excluding the EFO's assigned to the scattered 
positron. Bjorken-x (xgs) and W (Wes) are then defined by a combination of Ql and ?/e (eE 
method). 

The fractional momentum of the D*^ in the 7*p system is defined as 

.,^.) . M. P) 

where p*{D*) is the D*^ momentum in the 7*p center-of-mass frame. For the boost to the 7*p 



system, the virtual-photon vector was reconstructed using the double-angle (DA) estimator 
of the scattered positron energy, E'j^^. In this method, only the angles 9'^ and are used |]25[ . 
E'^^ is less sensitive to radiative effects at the leptonic vertex than the scattered e"*" energy 
determined using the electron method. 

For the KAn analysis, Q^, x and W were determined using the DA estimators {Qj^j^yXoA 
and Wda)- 

5 Monte Carlo simulation 



A GEANT 3.13-based ||26| Monte Carlo (MC) simulation program which incorporates the best 
current knowledge of the ZEUS detector and trigger was used to correct the data for detector 
and acceptance effects. The event generator used for the simulation of the QED radiation 
from the leptonic vertex was RAPGAP ^ interfaced to HERACLES 4.1 The charm 

quarks were produced in the BGF process calculated at leading order (LO). The charm mass 
was set to 1.5 GeV. The GRV94HO parton distribution functions (pdf's) were used for the 
proton. Fragmentation was carried out using the Lund model [EO] , as implemented in JETSET 



7.4 ISlj], with the full parton shower option. The fraction of the original c-quark momentum 



which is carried by the D*^ is determined from the 'SLAC fragmentation function, which is 



equivalent to the Peterson model |^2|], with the fragmentation parameter e set to 0.035 ||33 
The HERWIG 5.9 event generator was also used, with the same pdf's and c-quark mass as 
used in RAPGAP, to investigate the effects of fragmentation. 

Generated events with at least one D*^ — > D^nf [K^-k^ or K^'k^'k^'k^)'k'^ (or c.c.) 
were selected. These events were then processed through the detector and trigger simulation 
and through the same reconstruction program as was used for the data. 

6 Event selection 
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6.1 Trigger 



Events were selected online with a three-level trigger |T^. At the first level (FLT), inclusive 



DIS events are triggered by the presence of an isolated electromagnetic cluster in the RCAL 
or any energy deposition in excess of 3 GeV in any electromagnetic section of the CAL ||35|| . 
During high-luminosity periods, when the rate was high, a coincidence with an FLT track was 
also required. Tracks at the FLT are defined as a series of CTD hits pointing to the nominal 
interaction point. The efficiency of this trigger, with respect to the calorimeter-only trigger, 
was greater than 99.5% and in good agreement with the MC simulation. 

At the second level, algorithms are applied to reduce the non-e'^p background. The full 
event information is available at the third level trigger (TLT). At this level, events are accepted 
as DIS candidates if a high-energy scattered positron candidate is found within the CAL ('in- 
clusive DIS trigger'). Because of the high rate of low-Q^ events, this trigger was turned off in 
the region around the RCAL beampipe during high-luminosity operation. For the K27[ decay 
channel, a D*-finder (based on computing the Knn mass using tracking information and se- 
lecting loosely around the D*^ mass) was available in the TLT. Events at low were then 
kept by requiring a coincidence of an identified scattered positron anywhere in the CAL and a 
tagged Z)*^ candidate. This will be referred to as the 'D* trigger'. 

Using data selected from periods when both triggers were in use, the relative efficiency of the 
D* trigger with respect to the inclusive DIS trigger is found to be about 80% and independent 
of Q^, X, pt{D*) and ti{D*) within the measured kinematic regions. The MC simulations 
reproduce this efficiency to an accuracy better than the statistical accuracy of the data {'^ 
2%). 



6.2 Offline selection 

The DIS event selection was similar to that described in an earlier publication |^ ; namely, the 
selection required: 

• a positron, as identified by a neural network algorithm, with a corrected energy above 10 
GeV; 

• the impact point of the scattered positron on the RCAL was required to lie outside the 
region 26x14 cm^ centered on the RCAL beamline; 

• 40 < 5 < 65 GeV; and 

• a Z- vertex position {Z^t^l < 50 cm. 

The DIS events were restricted to the kinematic region 

• 1 < g2 < 600 GeV^ and 0.02 <y< 0.7. 

D*^ candidates were reconstructed from CTD tracks which were assigned to the recon- 
structed event vertex. Only tracks with at least one hit in the third superlayer of the CTD 
were considered. This corresponds to an implicit requirement that pt > 0.075 GeV. Tracks 
were also required to have |?7| < 1.75, where the pseudorapidity is defined as ?7 = — ln(tan|). 
The selected tracks were in the region where the CTD performance is well understood. For 
these tracks the reconstruction efficiency is above 95%. 

D*^ production was measured in the decay channel D*^ — > D^tt^ (+ c.c), which has a 
branching ratio of 0.683 ± 0.014 |^, followed by the decays — > K~7i^ (+ c.c.) or 
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K TT^TT TT^ 



[+ C.C.). The branching ratio for D° to Kn (KSn) is 0.0385 ± 0.0009 (0.076 
± 0.004) |3^. The remaining selection criteria were different for the two final states and are 
discussed separately. 



6.3 Selection for the K27r final state 

Pairs of oppositely-charged tracks were first combined to form a candidate. Since no particle 
identification was performed, the tracks were alternatively assigned the masses of a charged kaon 
and a charged pion. An additional slow track, with charge opposite to that of the kaon track 
and assigned the pion mass (tt^), was combined with the candidate to form a D*^ candidate. 
The combinatorial background for the K27t decay channel was further reduced by requiring 

• that the transverse momenta of the K and the vr were greater than 0.4 GeV, and that of 
the TTg was greater than 0.12 GeV. 

In addition, the momentum ratio requirement 

• p(D0)/j9(7r,) > 8.0 

was imposed. This requirement was used in the D* trigger to reduce the rate of candidate 
events with large AM = {Mk-k-k^ — Mk-k)-, far from the signal region. 
The D*^ kinematic region of the present analysis was defined as 

• 1.5 < Pt{D*) < 15 GeV and \ri{D*)\ < 1.5. 

Finally, the signal regions for the mass of the candidate, M(Z)°), and AM were 

• 1.80 < M(D°) < 1.92 GeV, and 

• 143 < AM < 148 MeV. 



6.4 Selection for the K47r final state 

Permutations of two negatively- and two positively-charged tracks were first combined to form 
a D'^ candidate. As for the i^27r channel, the tracks were alternatively assigned the masses 
of a charged kaon and a charged pion. An additional track, with charge opposite to that of 
the kaon track and assigned the pion mass (vr^), was combined with the candidate to form 
a -D*^ candidate. The combinatorial background for the KAtt decay channel was reduced by 
requiring 

• that the transverse momentum of the K and each tt was greater than 0.5 and 0.2 GeV, 
respectively, and that of the Hg was greater than 0.15 GeV. 

In addition, the momentum ratio requirement 

• p{D')/p{ns) > 9.5 

was imposed. 

The D*^ kinematic region was defined as 

• 2.5 < priD*) < 15 GeV and \viD*)\ < 1.5. 

The signal regions for M(D°) and AM = (MKmrmr^ — Mk-k-k-k) were 

• 1.81 < M(L)O) < 1.91 GeV, and 

• 143 < AM < 148 MeV. 
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6.5 Mass distributions 



Figures H(a) and (c) show the distributions of M{D^) for candidates with AM in the signal 
region for the two final states, while Figs, ^(b) and (d) show the distributions of AM for 
candidates with M{D^) in the signal region. Clear signals are observed around the expected 
mass values. 

For the K2tt final state, a fit to the M{D^) distribution of two Gaussians plus an exponen- 
tially falling background gives a peak at M(L>°) = 1863.2 ±0.8 MeV and a width of 23 ±2 MeV. 
The second Gaussian around 1.6 GeV in Fig. |^(a) originates primarily from decays to 
K~TT'^TT^ in which the neutral pion is not reconstructed. For the KAir channel, the background 
level was determined by using the side-bands outside the AM signal region (see the dashed 
histogram in Fig. 0(c)) to make a M{D^) distribution. The fit to the M{D^) distribution, made 
by adding a Gaussian to the background distribution, yielded M{D^) = 1862.7 ± 1.5 MeV and 
a width of 20 ± 2 MeV. The deviations of the data from the fit, visible in the region just above 
the signal in Fig. H(c), are mostly due to the mass misassignments of the K and vr candidates 
with the same charge from decay. This was verified by MC studies |3^. The mass values 
found for the are consistent with the PDG ||3^ value of 1864. 6±0. 5 MeV. 

The solid curve in Fig. 0(b) shows a binned maximum-likelihood fit to the AM distri- 
bution from the K27i channel using a Gaussian plus a background of the form A{AM — 
m^)^exp[C(AM — m,r)], where A, B and C are free parameters and is the pion mass. 
The fit to the K2tt plot gives a peak at AM = 145.44 ± 0.05 MeV, in good agreement with 
the PDG value of 145.397 ± 0.030 MeV, and a width of 0.79 ± 0.05 MeV, in agreement with 
the experimental resolution. The multiplicative exponential term is needed to describe the 
background suppression at large AM, which comes from the requirement on the momentum 
ratio p{D'^)/p{'7is). 

The solid curve in Fig. |](d) shows a fit of the AM distribution from the KAn channel using 
a Gaussian plus a background distribution obtained from the side-bands outside the M{D^) 
signal region (see the dashed histogram in Fig. §(d)). The fit yields a peak at AM = 145. 61± 
0.05 MeV and a width of 0.78 ± 0.07 MeV. 

For the K27r channel, the number of D*^ events obtained from a fit to the AM distribution 
in the restricted region of Q^, y, pt{D*) and f]{D*) is 2064 ± 72. The number of events in the 



KAti channel is determined from the AM distribution using the side band method [14|, which 



properly accounts for the combinatorial background and the background arising from the mass 
misassignments. The side bands, 1.74 < M{Ktx'K'k) < 1.79 GeV and 1.93 < M{K'K'K'k) < 
1.98 GeV, were normalized to the AM distribution in the region 150 < AM < 160 MeV. The 
number of D*^ events in the K4tt channel is found to be 1277 ± 124 in the restricted kinematic 
region of the measurement. 



7 Data characteristics 

The properties of the selected events are compared with those of the RAPGAP Monte Carlo 
simulation. All distributions shown are background-subtracted since they represent the number 
of signal events obtained by fitting the various mass distributions in a given bin. The data for 
the K27r channel are shown in Figs. |^ and ^ as solid points and the RAPGAP simulation as 
shaded histograms. All MC plots are normalized to have the same area as the data distributions. 

Figure |^ shows histograms of E'^, 9'^, jh and S and Fig. ||(a-c) displays the distributions of 
Ql, XeT. and VFeS- In general, reasonable agreement is observed between data and the MC 
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simulation. Figure ||(d-f) shows the transverse momentum, pt{D*), the pseudorapidity, ri{D*), 
and the energy fraction carried by the D*^ in the 7*p center-of-mass frame, x{D*). Ahhough 
the pt{D*) spectrum of the data is well described, the MC pseudorapidity spectrum is shifted 
to lower rj compared to the data and the x{D*) spectrum for the MC is shifted to slightly larger 
values. These discrepancies are examined in more detail below. 

The HERWIG [0] Monte Carlo was used for systematic studies. This MC describes the 
D*^ differential distributions better than RAPGAP, but does not give as good agreement with 
the DIS variables as RAPGAP since it does not contain QED radiative effects. 

Photoproduction, where the final positron is scattered through very small angles and escapes 
undetected through the RCAL beamhole, is a possible background source. Hadronic activity in 
the RCAL can be wrongly identified as the scattered positron, giving rise to fake DIS events. 
The effect was investigated using a large sample of photoproduced Z)*^ events generated with 
the HERWIG MC. After the final selection cuts the photoproduction contamination is found to 
be less than 1%, much smaller than the statistical error of the measurement, and so is neglected. 

The overall contribution to D*^ production from b quark decays in the measured kinematic 
region is estimated to be less than 2%, using HVQDIS with rrih = 4.5 GeV and a hadronization 
fraction oih ^ D* of 0.173 |3^. A similar study using RAPGAP yields an estimate of ~1% 



at low Q and less than 3% at high Q . Hence the contribution from h quark decays has been 
neglected. 



8 Systematic uncertainties 

The experimental systematic uncertainties in the cross section are grouped into several major 
categories: 

• Systematic uncertainties related to the inclusive DIS selection of the events: variations 
were made in the y cut, the RCAL box cut, and the vertex position cut. In addition, for 
the K2t: final state both and y were determined using the DA method rather than 
using the eS method. The combined variations resulted in a change of ±1.5% to the 
nominal cross section. For the KAtt channel, the electron method was used rather than 
the DA method and the combined variations resulted in a change of ±4.7%. 

• Systematic uncertainties in the D*^ selection: for the K2Ti final state, the minimum 
transverse momentum of tracks used in the D*^ reconstruction was raised and lowered 
by up to 100 MeV for the K and vr and by 25 MeV for the tt^ and resulted in a ±4.5% 
variation. The momentum ratio p{D^)/p{tTs) was raised by ±0.5 and yielded a negligible 
change. For the KAir channel, similar changes in the minimum transverse momenta and 
varying p{D^)/p{tTs) by ±0.5 combine to give a ±8.5% variation. 

• Systematic uncertainties related to the estimation of the number of events and back- 
ground uncertainty: for the K27r analysis, a AM-fit instead of a (binned) logarithmic- 
likelihood fit was performed. The M(D*') signal range, within which the AM distributions 
were fitted, was varied by ±10 MeV. These variations resulted in a change of ±2.5%. For 
the K4n analysis, the M[D^) signal range was also varied by ±10 MeV and the width of 
the side-bands used to estimate combinatorial background was varied. These variations 
resulted in a change of ±5.8%. 

• The systematic uncertainty related to the MC generator was estimated for the K2n anal- 
ysis by using the HERWIG MC generator to calculate the corrections for the cross section 
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determination. This variation yields a change of -1.3%. The larger overall systematic un- 
certainty for the KAn channel means that this systematic uncertainty was negligible for 
this channel. 



Since approximately 10% of the D*^ events |^0| are produced through a diffractive mech- 



anism which was not included in the Monte Carlo generators used to correct the data, 
acceptance corrections have also been obtained using a sample of diffractive events gen- 
erated with RAP GAP. The difference in the global correction factor for the diffractive 
events was less than 10%. This yielded a 1.3% variation in the overall cross section and 
was neglected. 

• The systematic uncertainty related to the trigger for the K27t final state was estimated 
to be ±1% by an analysis using only the inclusive DIS triggers and was neglected. 

• The overall normalization uncertainties due to the luminosity measurement error of 
±1.65%, and those due to the D*^ and decay branching ratios were not included 
in the systematic uncertainties. 

The systematic uncertainties were added in quadrature. The total systematic uncertainty is 
less than the statistical error for most of the differential distributions and is of the same order 
as the statistical error for the integrated cross section. 



9 Cross sections 

The cross sections for a given observable Y were determined from the equation: 



dY A-C-B-AY ' ' 

where A^ is the number of D*^ events in a bin of Ay, A is the acceptance (including migrations, 
efficiencies and radiative effects) for that bin, C is the integrated luminosity and B is the product 
of the appropriate branching ratios for the D*^ and D^. 

The RAPGAP MC was used to estimate the acceptance. In the K27r (KAir) kinematic region 
the overall acceptance was 25.4 (18.0)%. The statistical error of the MC is negligible compared 
to that of the data. 

The measured cross section in the region 1 < < 600 GeV^, 0.02 < y < 0.7, for 1.5 < 
Pt{D*) < 15 GeV and |?7(-D*)| < 1.5 using the K2n final state, is 

a{e+p e+D*^X) = 8.31 ± 0.31(stat.) j:°JO(syst.) nb, 

and for 2.5 < pt{D*) < 15 GeV and \i]{D*)\ < 1.5 the cross section using the KAtt channel is 

a{e+p e+D*^X) = 3.65 ± 0.36(stat.)l|]i?(syst.) nb. 

These results are in good agreement with the HVQDIS calculations of 8.44 and 4.13 nb, 
respectively. The parton distribution functions resulting from a ZEUS NLO QCD fit to 
the ZEUS, NMC and BCDMS data were used in these calculations. In this fit, only the gluon 
and three light quark flavors were assumed to be present. HVQDIS fragments the charm quark 
to a D*^ using the Peterson fragmentation function. For the above calculation, rric was set 
to 1.4 GeV, the Peterson fragmentation parameter e=0.035, and the mass factorization and 
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renormalization scales were both set to \J^m1 + Q^. The hadronization fraction /(c D*^) 
was set to 0.222 ± 0.014 ± 0.014, as measured by the OPAL collaboration H. The result 
is, however, sensitive to the choice of the parameters. For example, varying rric from 1.3 to 
1.5 GeV results in a variation of ±0.55 nb to the K2'n' prediction of 8.44 nb. If, instead of using 



the ZEUS NLO parton distributions, the three-flavor GRV98HO set is used, the calculated 
cross section for the i^27r kinematic region, assuming a charm mass of 1.4 GeV, is reduced from 
8.44 nb to 7.00 nb. 

Figures |^ and ^ show the D*^ differential cross sections in the measured kinematic regions. 
The data points are drawn at the bin center of gravity, which is defined as the point at which 
the value of the assumed theoretical curve (HVQDIS with central choice of parameters and 
RAPGAP fragmentation, see the discussion below) equals the mean value of the curve in the 
bin. The cross sections are compiled in Table ^ The 9% error on f{c—>- D*^) introduces an 
extra normalization uncertainty on the theoretical predictions which is not shown. 

The data are compared with two different theoretical calculations. There is good agreement 
with HVQDIS, shown as the open bands in Fig. |], except for the rj{D*) and x{D*) distribu- 
tions, where the HVQDIS calculations show a shift with respect to the data to more negative 
f]{D*) values and larger x{D*) values. 

The disagreement between the ri{D*) distribution and HVQDIS has also been observed by 
the HI collaboration in their DIS charm production data 0|. A large discrepancy in the 



Ti{D*) distribution is also observed between the ZEUS photoproduction of charm and a 



massive NLO calculation The discrepancy in shapes between the data and the HVQDIS 
prediction could result from the use of the Peterson fragmentation function. This predicts 
the magnitude of the momentum of the D*~^ from c — > D*~^ but produces no D*'^ transverse 
momentum relative to the c quark. In addition, no QCD evolution is included. In contrast, 
the fragmentation models in JETSET and HERWIG predict a migration of the charm quark 
emerging from the hard interaction towards positive pseudorapidities as it fragments into a 
D*~^ due to the interaction between the color charges of the c quark and the proton remnant. 



This has been called the 'beam-drag effect' |44 



To quantify the contribution of these fragmentation effects to the predicted D*^ differential 
cross sections, the RAPGAP Monte Carlo charm quark distribution was reweighted to match 
the two dimensional Pt{c), ri{c) distributions for the charm quarks from HVQDIS. The frag- 
mentation of these quarks was simulated by the Monte Carlo program including parton shower 
effects^. This reweighting procedure leads to an improved description of the data, as can be 
seen in Fig. ^ (shaded bands), where the two different implementations of the c — >■ D*~^ frag- 
mentation are compared. In order to study systematic effects, the same procedure as above was 
followed but with the parton shower option switched off and HERWIG used. This produced 
qualitatively similar results to the nominal method^. The data are again compared to the 
HVQDIS calculation including the RAPGAP-based fragmentation in Fig. ^. The solid curves 
are from using rric = 1.4 GeV, corresponding to the central mass value of the shaded band in 
Fig. i 

The data are also compared in Fig. | with the BKL tree-level calculations (dashed 
curves) in which the D*^ is created from both color singlet and color octet contributions to the 

^Note that this procedure keeps the cross section in the overall phase space in pt{D*) and rj{D*) equal to 
the HVQDIS result and is equivalent to applying the RAPGAP fragmentation to the HVQDIS charm quark 
prediction. 

Some double counting may occur between the NLO calculations and the parton shower, which contains 
resummcd terms of all orders. This effect is estimated to be small from the calculations with the parton shower 
option in JETSET switched off. 
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(cg)-state. The relative weight of these color configurations has been taken from comparisons 
of the BKL calculation to the published ZEUS D* cross sections in photoproduction The 
calculation shows reasonable agreement with these DIS data. 

Tables ^ and ^ give the resulting integrated cross sections binned in and y for the K2t: 
and A'47r final states, respectively. The bin widths were chosen such that they contained of the 
order of 100 signal events. The resulting purity in each bin is better than 70%. The resolutions 
in both and y are better than 10% in all bins. These bins were chosen to measure F^'^ 
as described in the next section. Tables § and ^ also show the HVQDIS predictions for the 
different kinematic bins. The predictions are in good agreement with the data. 

The quantitative agreement obtained between data and the HVQDIS calculation displayed in 
Fig. ^ and Tables ^ and ^ represents a confirmation of the hard scattering factorization theorem, 
in that the same gluon and three light-quark-fiavor parton distributions describe both the ZEUS 
F2 data and the D*^ differential cross sections reported here. In view of the agreement observed 
here, the HVQDIS program can be used to extrapolate outside the accessible kinematic region 
to obtain the total D*^ cross section. 



10 Extraction of 

The charm contribution, ¥2'^, to the proton structure function F2 can be related to the double 
differential cc cross section in x and by 
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.{[1 + (1 - yfWix, Q') - y'Fl'ix, Q')}. (4) 



dxdQ^ xQ^ 

In this paper, the cc cross section is obtained by measuring the D*"^ production cross section 
and employing the hadronization fraction /(c D*^) to derive the total charm cross section. 
Since only a limited kinematic region is accessible for the measurement of D*^, a prescription for 
extrapolating to the full kinematic phase space is needed. The contribution of Fl^{x, Q'^) to the 
cross section in the measured Q^, y region is estimated from the NLO theoretical prediction 



to be less than 1% and is therefore neglected. Equation (4) defines F^^^ as arising from events 
with one or more charm particles in the final state, but it is not a unique theoretical definition. 
It depends on the scheme and parton distributions |]46 |. 



In order to measure the contribution of charm to the inclusive F2, the integrated cross sections 
in the and y kinematic bins of Tables ^ and ^ were extrapolated to the full pt{D*) and ti{D*) 
phase space using HVQDIS with the RAPGAP-based fragmentation corrections discussed in 
the previous section. Typical extrapolation factors for the K2-K [KAti) final state were between 
4 (10), at low Q^, and 1.5 (4), at high . This procedure neglects the possibility of additional 
contributions outside the measured region due, for example, to intrinsic charm . 



The extrapolated cross sections are converted into cc cross sections using the hadronization 
fraction of charm to D*+: f{c D*+) = 0.222 ± 0.014 ± 0.014 The use of this value 

from OPAL implicitly assumes that charm production in DIS and e"'"e~ annihilation produces 
the same fractions of the various charm final-states. The production of charm bound states, 
such as e~^p —>■ J/ipX, is not accounted for when using the LEP c D*^ branching fraction. 
However, the inelastic J/tp cross section has been calculated to be only 2.5-4.5% of the 



total charm production cross section predicted by HVQDIS in the Q^, y range of this analysis. 
The elastic J/ip cross section has been measured in DIS and is less than 0.5% of the 



predicted total charm cross section in the range of that measurement. The 9% uncertainty on 
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the c — * D*^ hadronization is larger than that arising from these J/ip contributions, which 
have consequently been neglected. 

The systematic uncertainty on the extrapolation of the measured D*^ cross sections to the 
full pt{D*) and r]{D*) phase space was investigated: varying the parameter rric by ±0.15 GeV 
gave a variation which was typically <5%; using the standard Peterson fragmentation parameter 
e = 0.035 (instead of the RAPGAP fragmentation correction) yielded changes typically <15%; 
and using the GRV98HO pdf's in the NLO calculation generally caused changes of <20%. 
If these uncertainties are added in quadrature, they are typically smaller than the statistical 
errors. However, the fact that the data are measured in a small part of the available phase space 
means that a realistic uncertainty on the extrapolation cannot be evaluated. Therefore these 
extrapolation uncertainties are not included in the systematic uncertainties discussed below. 

Since the structure function varies only slowly, it is assumed to be constant within a given 
and y bin, so that the measured F^^^ in a bin i is given by 

<^i,theor\^ V U y^) 

where the cross sections CTj in bin i are those for the measured pt{D*) and f]{D*) region, and 
the subscripts meas and theor denote 'measured' and 'theoretical', respectively. The value of 
^Ttheor calculated from the NLO coefficient functions |]^, as implemented in a convenient 



parametrization ||4^. The functional form of F2'^theor used to quote the results for F2'^ at 
convenient values of Xi and close to the center-of-gravity of the bin. In this calculation, 
the same parton densities, charm mass (mc= 1.4 GeV), and factorization and renormalization 

scales {^Aml + Q"^) have been used as for the HVQDIS calculation of the differential cross 
sections. 



10.1 Combination of F^^"^ from both decays 

Finally, the results from the two decay channels were combined in the eight common bins, tak- 
ing into account all systematic uncertainties. The combined value is a weighted average with 
weights according to the statistical precision of the individual measurements. The systematic 
uncertainties were assumed to be either uncorrelated or 100% correlated between the analy- 
ses, as appropriate. All uncertainties concerning the DIS event selection were assumed to be 
correlated. Only the effect of the variation of the mass window and the changes in the pt 
requirements for the decay products were taken as uncorrelated. As in both analyses, the 
positive and negative errors were treated separately. The procedure leads to a gain in statistical 
precision of 5-25%, compared to using only the K2tt decay channel. 



10.2 Results and discussion 

Table ^ and Fig. ^ display the Fl'^ values in the various bins as a function of x. The structure 
function shows a rise with decreasing x at constant values of . The rise becomes steeper 
at higher Q^. 

The curves in Fig. ^ represent the results of the NLO QCD calculation |^ with the ZEUS 



NLO QCD pdf's. The central, solid curve corresponds to a charm quark mass of 1.4 GeV. 
Since good agreement was obtained between data and the HVQDIS calculation for the D*"^ 
differential cross sections and for the integrated cross sections shown in Tables |^ and ^ and 
since that calculation was used to extrapolate to the full kinematic range, the curves would be 
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expected to describe the resulting values of F^'^. The total uncertainty in the calculation of 
F^, shown as the band of dashed curves around the solid curve, corresponds to the uncertainty 
propagated from the ZEUS NLO QCD fit and is dominated by the uncertainty in the charm 
quark mass, which was varied from 1.2 to 1.6 GeV. 

Figure |^ shows Fg'^ at constant x values as a function of . Although the number of points 
is small, large scaling violations of the structure function are evident. The curves superimposed 
on the data are from the same calculation as shown in Fig. |^. 

Figure ^ shows the ratio of ^2^^ to F2, the inclusive proton structure function, as a function 
of X in fixed-Q^ bins. The curves superimposed on the data are again from the calculation used 
for Fig. 1^. The values of F2 used to determine the ratio were taken from the ZEUS NLO QCD 
fit at the same and x for which F!^"^ is quoted. The error on F2 is negligible in comparison 
to F2^. The charm contribution to F2 rises steeply with decreasing x. In the measured x 
region, F2'' accounts for < 10% of F2 at low and x ~ 5 • 10""^ and rises to ~ 30% of F2 for 
11 GeV^ at the lowest x measured. The strong rise of ^2^*^ at low values of x is similar to 
that of the gluon density and thus supports the hypothesis that charm production is dominated 
by the boson-gluon-fusion mechanism. 

11 Summary 

This paper presents an analysis of D*^ production in DIS using the combined ZEUS 1996 and 
1997 data samples with an integrated luminosity of 37 pb~^, about ten times larger than in 
the previous ZEUS study. In addition, both the K2'k and KAtt decay modes of the D* have 
been employed and their results combined. In the experimentally accessible region of 1.5 (2.5) 
<PT{D*)<lh GeV and \'r]{D*)\ < 1.5, the cross section for D*^ production for the K2'k {KA-k) 
final state in events with 1 <Q'^< 600 GeV^ and 0.02 < ?/ < 0.7 is 8.31 ±0.31(stat)t[J|[j(sys) nb 
(3.65 ± 0.36(stat)l[5:|?(sys) nb). 

QCD calculations of charm production based on the NLO boson-gluon-fusion process with 
three flavors of light quarks show excellent agreement with the overall cross section and with the 

and y distributions. The ti{D*) and x{D*) distributions, however, cannot be reproduced with 
the standard Peterson fragmentation. Good agreement is obtained after a more appropriate 
c D*^ fragmentation, such as that in JETSET, is used. 

The quantitative agreement between the NLO pQCD calculations and the ZEUS data pro- 
vides a confirmation of the hard scattering factorization theorem, whereby the same gluon 
density in the proton describes both the inclusive F2 and the DIS production of charm. 

The charm contribution, ^2^^^, to the proton structure function F2 was obtained using the NLO 
QCD calculation to extrapolate outside the measured pt{D*) and f]{D*) region. Compared to 
the previous ZEUS study, the kinematic range has been extended down to = 1.8 GeV^ and 
up to = 130 GeV^, with reduced uncertainties. The structure function F2'^ exhibits large 
scaling violations, as well as a steep rise with decreasing x at constant Q^. For Q'^> 11 GeV^ 
and X ~ 10~^, the ratio of F2'^ to F2 is about 0.3. 
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Table 1: The D*^ KrcrCs and D*^ KTnnrirs differential cross sections. The bin range, the 
center- of- gravity of the bin (see text) and the cross sections for all the data in Figs. and^ are 
shown. The first error is the statistical error and the asymmetric errors are the statistical and 
systematic uncertainties added in quadrature. The overall normalization uncertainties arising 
from the luminosity measurement (±1.65%) and from the D*^ and decay branching ratios 
are not included. 



KnTTg 


KnTTg 


(range) logioiQ"^) 


da / dlogiQ^Q"^) 
(nb) 


(range) logio{x) 


da / dlogio{x) 
(nb) 


(0.0,0.7) 0.39 
(0.7,1.0) 0.85 
(1.0,1.3) 1.16 
(1.3,1.6) 1.45 
(1.6,1.9) 1.74 
(1.9,2.3) 2.08 
(2.3,2.8) 2.48 


5.99±0.47 tlfi 
5.17±0.39 
4.50±0.33 tlfi 
2.72±0.25 
1.47±0.17 
0.47±0.08 
0.135±0.073 1[J:|]?^ 


(-4.1,-3.4) -3.69 
(-3.4,-2.8) -3.08 
(-2.8,-2.3) -2.56 
(-2.3,-2.0) -2.16 
(-2.0,-1.5) -1.82 


4.54 ± 0.35 tofe 
4.24 ± 0.23 toir 
2.06 ± 0.16 tm 
0.78 ± 0.14 tin 
0.25 ± 0.10 tin 


KtttTs 




(range) W 
(GeV) 


da/dW 
(nb/GeV) 


(range) x{D*) 


da/dx{D*) 
(nb) 


(50,90) 73 
(90,115) 102 
(115,145) 129 
(145,175) 159 
(175,200) 187 
(200,250) 222 


0.0450±0. 0035+^:^^1^ 
0.0659±0.0048l[J:|]|]^^ 
0.0510±0.0039i[;:l][Ji 
0.0442±0.0038l[5:|]|ji^ 
0.0361±0.0043l[5:|]|]?I 
0.0266±0.0032t[J:|]|]i 


(0.0,0.2) 0.13 
(0.2,0.35) 0.28 
(0.35,0.5) 0.42 
(0.5,0.6) 0.55 
(0.6,0.75) 0.67 
(0.75,1.0) 0.80 


9.3 ± 1.2 t2*2 
13.2 ± 1.2 t\i 
13.61 ± 0.98 tlfo 
12.79 ± 1.03 tin 
8.15 ± 0.60 t^iil 
1.12 ± 0.10 tin 






(range) pt{D*) 
(GeV) 


da/dpT{D*) 
(nb/GeV) 


(range) Pt{D*) 
(GeV) 


da/dpT{D*) 
(nb/GeV) 


(1.5,2.4) 1.91 
(2.4,3.1) 2.72 
(3.1,4.0) 3.50 
(4.0,6.0) 4.77 
(6.0,15.0) 7.93 


3.82±0.39 tlf^ 
2.72±0.21 toll 
1.57±0.10 till 
0.527±0.033 ttlilfi 
0.0419±0.0034j:[J:J][J^? 


(2.5,3.0) 2.74 
(3.0,3.5) 3.24 
(3.5,4.0) 3.73 
(4.0,6.0) 4.74 
(6.0,8.0) 6.77 
(8.0,10.0) 8.75 
(10.0,15.0) 11.73 


2.42 ± 0.76 tlfr 
1.63 ± 0.49 tlfi 
1.03 ± 0.29 tin 
0.43 ± 0.06 tin 

0.127 ± o.oi9;j]:[!^^ 
0.044 ± o.oiotJ]:[!JJ 
0.007 ± o.oo2lJ]:[![!^ 


KtttTs 




(range) rj{D*) 


da/dT]{D*) 
(nb) 


(range) rj{D*) 


da/dr]{D*) 
(nb) 


(-1.5,-0.8) -1.13 
(-0.8,-0.35) -0.58 
(-0.35,0.0) -0.18 
(0.0,0.4) 0.20 
(0.4,0.8) 0.60 
(0.8,1.5) 1.15 


1.93±0.19 tliil) 
2.41±0.21 
3.04±0.28 toij 
3.09±0.25 toil 
3.15±0.27 tin 

3.26±o.3i tin 


(-1.5,-1.0) -1.18 
(-1.0,-0.5) -0.73 
(-0.5,0.0) -0.24 
(0.0,0.5) 0.26 
(0.5,1.0) 0.75 
(1.0,1.5) 1.22 


0.73 ± 0.21 tlfi 
1.03 ± 0.23 tm 
1.48 ± 0.26 toll 
1.48 ± 0.31 tin 
1.12 ± 0.35 tin 
2.08 ± 0.45 tlfs 
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Table 2: The cross sections for L>* production from the K2'k final state. The table contains 
for each bin: the range of the bin; the y range of the bin; the measured D*^ cross section 
in the bin with statistical and systematic uncertainties; and the HVQDIS prediction for this 

cross section. The overall normalization uncertainties arising from the luminosity measurement 
(±1.65%) and from the D*^ and decay branching ratios are not included. 



range 

V ) 


y 

range 


rnpns ( nil 1 


a{D*) 
nred i nl) 1 


1-3.5 


0.70-0.24 
0.24-0.11 
0.11-0.02 


1.45 ±0.231^:^^ 
1.00±0.20i[]i? 

92 ± 16+nnQ 


1.12 
1.00 
1.10 


3.5-6.5 


0.70-0.22 
0.22-0.11 

0.11-0.02 


0.73 ± 0.10l^:^{ 
0.342 ± 0.055^^36 

433 ± 066+nni^« 


0.54 
0.405 

0.528 


6.5-9 


0.70-0.15 
15-0 02 


0.388 ± 0.060^^:^^^ 
288 ± 046"^nnni 


0.369 
0.355 


9-14 


0.70-0.23 
0.23-0.11 

0.11-0.02 


0.370 ± 0.0571^:^^^ 
0.314 ± 0.0451° 076 
0.253 ± 0.042l[]:[!^J 


0.302 
0.251 

0.316 


14-22 


0.70-0.23 
0.23-0.11 
0.11-0.02 


0.25 ± 0.051^:^^ 
0.254 ± 0.043lH^^ 
0.226 ± 0.035lg:|]^^ 


0.26 
0.212 
0.250 


22-44 


0.70-0.23 
0.23-0.11 
0.11-0.02 


0.387 ± 0.0531^:^^^ 
0.200 ± 0.027lg:[]^? 
0.198 ± 0.033l|j:|]?^ 


0.301 
0.226 
0.240 


44-90 


0.70-0.23 
0.23-0.02 


0.202 ± 0.0431^:^?^ 
0.200 ± 0.032l[|:[]l^ 


0.188 
0.221 


90-200 


0.70-0.23 
0.23-0.02 


0.090 ± 0.0231^:^^^ 
0.075 ± 0.015tHo7 


0.099 
0.086 
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Table 3: The cross sections for D*^ production from the KAn final state. The table contains 
for each bin: the range of the bin; the y range of the bin; the measured D*^ cross section 
in the bin with statistical and systematic uncertainties; and the HVQDIS prediction for this 
cross section. The overall normalization uncertainties arising from the luminosity measurement 
(±1.65%) and from the D*^ and decay branching ratios are not included. 



range 
(GeV2) 


y 

range 


a{D*) 
meas. (nb) 


a{D*) 
pred. (nb) 


1-10 


0.70-0.34 
0.34-0.02 


0.94±0.29+^:^^ 
1.39±0.5lt[J:g^ 


0.61 
1.55 


10-21 


0.70-0.28 
0.28-0.02 


o.35±o.o8i;]:i(t 

0.59±0.1ltHI 


0.21 
0.47 


21-33 


0.70-0.22 

0.22-0.02 


0.115±0.055;n23 

0.125±0.045i[]:H 


0.154 

0.209 


50-600 


0.70-0.22 
0.22-0.02 


o.23±o.o7+^:y^ 

0.240±0.06ll[i[Ji 


0.249 
0.221 
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Table 4: The F2'' results derived from the combination of the K2tt and KAn channels (see 
teoct). The table contains for each bin: the value at which Fg*^ is reported; the x value at 
which ^2^^ is reported; and the measured Fg'^ with statistical and systematic uncertainties. The 

overall normalization uncertainties arising from the luminosity measurement (±1.65%), the 
D*^ and decay branching ratios, the charm hadronization fraction to D*^ (±9%) and the 
extrapolation uncertainties (see text) are not included. 
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Figure 2: 2'£^CAS' c/ato in the kinematic region 1 < < 600 GeV'^and 0.02 < y < 0.7; (a) 
M{D^) = M(X7r) anc? (h) AM = M(i^7r7r,) - M{Ktt) for the K2tt final state. The data are 
shown in the same kinematic region for the KAtt final state in (c) M{D^) = M{K3n) and (d) 
AM = MiKSTTTTs) - M{K?,tt). The events in (a) and (c) are those for which 1^3 < AM< I48 
MeV. Similarly, (b) and (d) are shown for events in the M{D^) signal region: 1.80-1.92 (1.81- 
1.91) GeV for the K2t: (Kin) channel. The solid curves are the results of the fits described in 
the text. The dashed histograms in (c) and (d) represent the background distributions obtained 
by restricting the sample to side-bands from the AM and M{D'^) distributions, respectively. 
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Figure 3: A comparison of reconstructed DIS quantities at the detector level for the K2ti data 
(points) and for the RAPGAP Monte Carlo simulation (shaded histogram): (a) the scattered 
positron energy, E'^, (b) the scattered positron angle, 9'^, (c) the hadronic angle, 7/j and (d) 
6= T.i{Ei -pz,i)- 
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p^(D*) (GeV) ri(D*) x(D*) 

Figure 4: A comparison of the reconstructed DIS kinematic distributions at the detector level for 
the K2iT data (points) and for the RAPGAP Monte Carlo simulation (shaded histogram): (a) 
Ql, (b) Bjorken-x, Xes and (c) the total hadronic center-of-mass energy, Wes- The remaining 
plots show a comparison of the reconstructed D*^ -related quantities at the detector level for 
data and for the RAPGAP Monte Garlo simulation: (d) pt{D*), (e) r]{D*) and (f) x{D*). 
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Figure 5: Differential cross sections for D*"^ production from the K2Tr final state (solid dots) in 
the Q^, y, pt{D*) and rj{D*) kinematic region as functions of (a) log^^Q"^ , (h) log^g^; i^^) ^> 

(d) pt{D*), (e) rj{D*) and (f) x{D*). The inner error bars show the statistical uncertainties 
while the outer ones show the statistical and systematic uncertainties summed in quadrature. 
The results from the KAn channel (open triangles) are also shown in the pT{D*)(d) and ri{D*) 

(e) plots. The data are compared with the NLO QCD calculation as implemented in HVQDIS 
using the ZEUS NLO pdf's. The open band corresponds to the standard Peterson fragmentation 
function with the parameter e = 0. 035. For the shaded band, the Peterson fragmentation was 
replaced by that extracted from RAPGAP (see the text for details). The boundaries of the bands 
correspond to charm mass variations between 1.3 (upper curve) and 1.5 GeV (lower curve). In 
(a) and (b), the open band is indistinguishable from the shaded band. 
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Figure 6: Differential cross sections for D*^ production from the K2tt final state (solid dots) in 
the accessible , y, pt{D*) and rj{D*) region as functions of (a) logi^Q"^ , (h) log;^oX, (c) W , 
(d) pt{D*), (e) rj{D*) and (f) x{D*). The inner error bars show the statistical uncertainties 
while the outer ones show the statistical and systematic uncertainties summed in quadrature. 
The solid curves show the results of the HVQDIS calculation with RAPGAP-based fragmentation 
and rric = 1-4 GeVwhile the dashed curves correspond to the BKL results (see text). 
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Figure 7: The measured ^2*^ at values between 1.8 and 130 GeV^as a function of x. The 
inner error bars show the statistical uncertainty and the outer ones show the statistical and 
experimental systematic uncertainties summed in quadrature. The curves correspond to the 
NLO QCD calculation ^ using the result of the ZEUS NLO QCD fit to F2 The 
solid curves correspond to the central values and the dashed curves give the uncertainty due 
to the parton distributions from the ZEUS NLO fit. The overall normalization uncertainties 
arising from the luminosity measurement (±1.65%), the D*^ and decay branching ratios, 
the charm hadronization fraction to D*^ (±9%) and the extrapolation uncertainties (see text) 
are not included. 
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Figure 8: The measured at x values between 0.00005 and 0.02 as a function of . The 
various values of x are indicated to the right of the data points. For clarity of presentation, the 
F^'^ values have been scaled by the number shown in parentheses next to the x value. The inner 
error bars show the statistical uncertainty and the outer ones show the statistical and systematic 
uncertainties summed in quadrature. The curves correspond to the NLO QCD calculation j^, 
using the result of the ZEUS NLO QCD fit to F2 The solid curves correspond to the central 
values and the dashed curves give the uncertainty due to the parton distributions from the ZEUS 
NLO fit. Details of this calculation are given in the text. The overall normalization uncertainties 
arising from the luminosity measurement (±1.65%), the D*^ and decay branching ratios, 
the charm hadronization fraction to D*^ (±9%) and the extrapolation uncertainties (see text) 
are not included. 
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Figure 9: The ratio of to F2 at values between 1.8 and 130 GeV^as a function of x. The 
ratio is calculated using the F2'^ values measured in this analysis and the ZEUS NLO QCD fit 
to F2 l^j . The inner error bars show the statistical uncertainty and the outer ones show the 
statistical and systematic uncertainties summed in quadrature. The curves correspond to the 
NLO QCD calculation /0, usmg the result of the ZEUS NLO QCD fit to F2 /g7|/. The solid 
curves correspond to the central values and the dashed curves give the uncertainty due to the 
parton distributions from the ZEUS NLO fit. Details of this calculation are given in the text. 
The overall normalization uncertainties arising from the luminosity measurement (±1.65%), 
the D*^ and decay branching ratios, the charm hadronization fraction to D*^ (±9%) and 
the extrapolation uncertainties (see text) are not included. 
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